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Abstract 
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data sample collected by ALEPH from 1991 to 1995. The selection is based on the 
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decays, and on the topological properties of signal events. The result derived from the 
selected sample is g cS = (3.26 ± 0.23(stat) ± 0.42(syst))%. 
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1 Introduction 



In this letter a measurement of the production rate of cc pairs from gluons in hadronic 
Z decays is described. The selection relies on tagging semileptonic decays of the c 
quarks from gluon splitting, and makes use of several discriminating variables related 
to the event topology. 

The rate of gluon splitting to cc pairs is defined as 

N(Z -» qqg, g -» cc) 
9cZ N(Z -> hadrons) ' 1 ' 

Measuring g c5 is an important test of perturbative QCD at the Z scale. The 
processes g — > cc and g — > bb are also significant backgrounds for several analyses 
involving heavy quarks. For example, one third of the total experimental uncertainty 
on i?b comes from these processes pQ. Furthermore, gluon splitting to heavy quarks is 
a background for Higgs boson searches 0. 

The theoretical treatment of the production of heavy quarks from gluons is 
described in [HI IU E] and the rate g cZ is predicted to be in the range 1.4% to 2.5%. 

Previous measurements have been performed using a D* tag [HE], a lepton tag [312] 
or event shape variables The lepton tag currently provides the most precise results. 



2 The ALEPH detector 

A detailed description of the ALEPH detector can be found in ^H] and of its 
performance in [TT] . A brief overview is given in this section. 

Charged particles are detected in the central part of the apparatus, consisting of a 
high resolution silicon strip vertex detector (VDET), a cylindrical drift chamber (ITC) 
and a large time projection chamber (TPC). The three tracking detectors are immersed 
in a 1.5 T axial magnetic field provided by a superconducting solenoid. They are 
surrounded by the calorimetric system, consisting of the electromagnetic calorimeter 
(ECAL), the hadron calorimeter (HCAL) and the muon chambers. 

The VDET [TJ] lies at the core of the tracking system. It is made of two layers, at 
average radii of 6.5 and 11.3 cm, each providing measurements in both the r<j) and rz 
projections, with a resolution of 12 /xm for r<p coordinates and varying between 12 and 
22 fim for z coordinates, depending on the track polar angle 9. The angular coverage 
is |cos 9\ < 0.85 for the inner layer and |cos 9\ < 0.69 for the outer layer. 

The ITC measures up to eight coordinates per track in the r<p projection, with a 
resolution of 150 fim. 

The TPC provides up to 21 three-dimensional coordinates per track, with 
resolutions in the r<p and rz projections of 180 /mi and 500 /im, respectively. The 
TPC also provides up to 338 measurements of the specific energy loss by ionization 
(dE/dx); this allows electrons to be distinguished from other charged particles by more 
than three standard deviations up to a momentum of 8 GeV/c. 

The ECAL is a sampling calorimeter covering the angular range |cos 9\ < 0.98, 
segmented in 0.9° x 0.9° projective towers, read out in three longitudinal stacks. The 



1 



nominal thickness of the calorimeter is 22 radiation lengths. The energy resolution for 
isolated electrons and photons is <Je/ E = 0.009 + 0.18/V^E, with E measured in GeV. 

The hadron calorimeter (HCAL) is composed of 23 layers of streamer tubes 
interleaved with iron slabs. The total iron thickness corresponds to about 7 interaction 
lengths at normal incidence. 

Electrons are identified by the characteristic longitudinal and transverse 
development of their associated showers in the ECAL. The dE/dx information from 
the TPC is used to enhance the hadron rejection power, while non-prompt electrons 
originating from photon conversions in the detector material are rejected on the basis 
of their kinematical and geometrical properties. 

Muons are identified by their penetration pattern in the HCAL; the additional three- 
dimensional coordinates measured in two double layers of external muon chambers help 
in resolving the remaining possible ambiguities. 

The lepton identification technique is described in detail in \M\ , with minimum 
momentum cuts of 2 GeV/c for electrons and 2.5 GeV/c for muons. 



3 Preselection 

The analysis is based on the LEP 1 data set, which consists of about 3.9 million Z — ► qq 
decays collected by ALEPH from 1991 to 1995. The analysis makes also use of 8.7 
million simulated Z — > qq events, 5.1 million Z —>■ bb events, 2.3 million Z — > cc events, 
and 1.8 million signal events each containing the g —>■ cc process. The generator is 
based on JETSET 7.3 [T3], and all events are passed through a detailed simulation of 
the detector based on GEANT3 ^Hj- Simulated events are reweighted to take into 
account the latest world average of g cS and g h ^. 

The method relies on the analysis of events clustered into three jets, where an 
electron or a muon is found in the least energetic jet (taken to be the gluon jet 
candidate). Additional discrimination between signal and background is obtained using 
variables related to the event topology. 

First, events with a small value of the thrust are selected (T < 0.94). Figure [T] 
shows the thrust distribution for data and simulated hadronic events, together with 
the distribution of signal events. For the events selected, particles are clustered into 
three jets, using the JADE algorithm ^7] ("E scheme"). The energies of the jets are 
recalculated by enforcing energy and momentum conservation, under the assumptions 
that jet directions are perfectly measured and that jets are massless, as 

E% - E cm - ■ ■ ; : ■ ; — , [2 

sin ip jk + sin ipij + sin ip ik 

where E cm is the centre-of-mass energy and ip^ is the angle between jets i and 
j. Non-planar events, for which Eq. El does not hold, are rejected by requiring 
■012 + ^23 + ^31 > 358°. The jet energies calculated in this way are used to order 
the three jets by decreasing energy. 

Lepton candidates are searched for, following the same procedure as in |13j . In the 
present analysis the cut on the distance of closest approach of the lepton track to the 
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Figure 1: Thrust distributions for data and simulation, normalized to the same area. 



beam axis (do) is tightened to do < 0.1 cm, to improve the rejection of leptons from K 
or 7i decay. Events with an identified lepton belonging to the third jet are retained for 
further analysis. 

This preselection yields 13363 events, out of which 5639 contain an electron 
candidate and 7724 a muon candidate. The composition of the two subsamples 
evaluated with simulated events is shown in Table Q the flavour content of the 
background is substantially different in the two cases. The larger fraction of light 
quark events for the muon subsample is related to the contamination of pions. The 
fraction of different particle types contributing to the samples of muon and electron 
candidates is shown in Table where prompt lepton indicates a lepton originating 
from the decay of a heavy flavour particle. 



Table 1: Composition of the preselected sample. 





e 




g -> cc 


17.0 ± 0.3 % 


12.8 ± 0.3 % 


g — > bb 


3.4 ± 0.2 % 


2.5 ±0.1 % 


Z — > uu, dd, ss 


10.9 ± 0.3 % 


27.7 ±0.4 % 


Z — > cc 


21.5 ±0.4 % 


20.0 ± 0.3 % 


Z ^bb 


47.2 ± 0.4 % 


37.0 ± 0.4 % 
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Table 2: Fraction of different particle types in the samples of muon and electron 
candidates, after preselection. 



Muon candidates 


prompt [i 


56% 


71 — > \W 


25% 


misidentified tt 


10% 


K — > iw 


5% 


misidentified K 


3% 


others 


1% 



Electron candidates 


prompt e 


69% 


7 — > e + e~ 


26% 


misidentified tt 


4% 


others 


1% 



4 The extraction of g C c 

The purity of the selected sample is inadequate to perform a measurement of g cE with 
a meaningful precision. Additional variables are needed to provide further separation 
between the signal and the Z — > qq background, as follows. 

• After the preselection cut mentioned in Section as clearly visible from the 
distributions in Fig. ^ the thrust T retains some discriminating power between 
the (multi-jet-like) signal events and the bulk of the Z — > qq background. 

• Each event is divided into two hemispheres by a plane perpendicular to the thrust 
axis. The confidence levels Bi and B 2 (hereafter called b-tag probabilities [TH] ) 
that the charged particle tracks of each hemisphere originate from the primary 
vertex provide a tag against Z — ► bb events and, to a lesser extent, Z — > cc events. 

• The b-tag probability be and the momentum Pi of the tagged lepton in the third 
jet also contribute to the aforementioned anti-b-tag capability. 

• Conversely, the projection -P^ iss of the missing momentum along the axis of 
the third jet is largest for c — > I decays, and is therefore discriminant against 
Z — ► uu, dd and ss events, for which the missing momentum direction is mostly 
random. This variable is statistically almost as powerful as the third jet mass, 
used for instance in [7J E]. However, the latter is not included in the present 
analysis because it is found to depend on details of the fragmentation model and 
to be inadequately reproduced by the simulation. 

• Finally, the relative discriminating power of the variables mentioned above varies 
with the third jet polar angle. The cosine | cos 9(\ of the tagged lepton polar angle 
is therefore added to better control this dependence. 

These variables are combined with an artificial neural network into single 
discriminants, N e and N^, for the electron and the muon samples separately. The 
neural network is trained with half the simulated Z — > qq sample and half the simulated 
signal sample. The remaining statistics, together with the Z — > cc and the Z — > bb 
simulated samples, are used to determine the selection efficiencies. 



4 




Figure 2: Distributions of the signal and of the background components, normalized 
to the same area, for P^iss- 

As an example, the discriminant power of -P^iss is shown in Fig. where the 
distributions of the signal and background components are compared (normalized to 
the same area). Figure H3 shows the agreement between data and simulation for the 
variables with highest discriminating power, at preselection level. 

The distributions of iV e and are shown in Fig. |3J The separation between signal 
and background is worse for the muon subsample, due to the high contamination of 
non-prompt muons in the preselected sample. 

The cuts on N e and are chosen to be 0.55 and 0.45 in order to minimize the total 
uncertainty on g c5 . This final selection cut yields 2258 events with an electron candidate 
and 3332 events with a muon candidate. The sample composition is estimated from 
the simulation, and is shown in Table El 

The value of the gluon splitting rate is extracted as 

_ / ~ (1 ~ 9bb) e Q - 9bb e B /„\ 
9cc — , {O) 

e c - £q 

where / is the fraction of events selected in the data. The selection efficiencies 6q, 
ec, €b for events with, respectively, no gluon splitting, gluon splitting in cc, and gluon 
splitting in bb are estimated from the simulation (Table HJ). 

In the calculation of g c5 , g h ^ is fixed to the latest world average value 
(2.54 ± 0.51) x 10~ 3 PHI- The results extracted from the two subsamples are 
g% = (3.32 ± 0.28)% and = (2.99 ± 0.38)%, where the errors only account for the 
statistical uncertainty on /. 



5 





o 

CO 

> 

> 

pq 



10 3 b 




_i i i i_ 







_l I l_ 



> 

$2000 

CO 

J— 1 

§1000 
> 

pq 



-20 



: A 




— ^v* r *f*^*^ , i 


V 



20 



P^ ss (GeV/c) 



Figure 3: Distribution of the variables with highest discriminating power, in data and 
simulation. 
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Figure 4: Output of the neural network for electrons and muons. In these plots, the 
signal is normalized according to the results obtained from Eq. El 

Table 3: Composition of the selected sample, after the cuts on N e and iV M . 





e 




g -> cc 


25 % 


18 % 


g -> bb 


4 % 


2 % 


Z — > uu, dd, ss 


18 % 


35 % 


Z — > cc 


28 % 


22 % 


Z -> bb 


25 % 


23 % 



Table 4: Fraction of events selected and selection efficiencies for the different event 
categories, with statistical errors. 





electrons (%) 


muons (%) 


/ 


0.0600 ±0.0013 


0.0885 ±0.0015 




0.0431 ± 0.0006 


0.0747 ±0.0009 


ec 


0.491 ±0.005 


0.473 ± 0.005 


£b 


0.81 ±0.06 


0.80 ±0.06 
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5 Systematic errors 

The sources of uncertainty on the selection efficiencies given in Table 0] are discussed in 
this section. The resulting systematic uncertainties on g cE are summarized in Table |H1 

i. Statistics of the simulation 

The statistical uncertainties on the selection efficiencies translate to 
A^ E = ±0.13% and Ag& = ±0.21%. 

ii. Gluon splitting to bb pairs 

The uncertainty on g> b g yields Ag* s — ±0.09% and Ag^. = ±0.10%. 

iii. Thrust cut 

The efficiency of the thrust cut is found to be slightly but significantly different 
in data and simulation. The effect is taken into account by reweighting the 
simulated efficiencies to the data efficiencies. If no reweighting is performed, a 
shift is observed in the result: Ag^ = ±0.07% and Ag^ = ±0.11%. This shift is 
taken as an error estimate, with full correlation between the two channels. 

iv. Heavy quark properties 

The lifetimes of the b hadrons determine the impact parameter distribution of the 
decay products,which is the basis of the "anti-b tag" method used in this analysis. 
Their experimental values are taken from [20!; the uncertainties translate to 
A^l = ±0.03% and Ag^ = ±0.03%. 

The semileptonic branching ratios of c and b hadrons [21] affect the flavour 
composition of the samples. Their uncertainties and their influence on the result 
of the analysis are shown in Table 

The energy spectra of leptons coming from b — ► £, c — > £ and b — > c — > £ decays 
have been tuned and varied as in [221, an d the effect has been propagated to the 
measured gluon splitting rate as shown in Table El 

The xe (= PhadZ-E-beam) distribution of the heavy hadrons affects their 
decay length distribution and the momenta of their decay products. The 
uncertainty coming from the limited knowledge of this quantity has been 
evaluated by varying the mean and x^ within their estimated errors: 

Table 5: Effect of the uncertainties on the rates and energy spectra of heavy hadron 
semileptonic decays. 





As c c E (%) 


A^ E (%) 


BR(b -> I) = (10.65 ± 0.23)% 
BR(b -f c -f I) = (8.04 ± 0.19)% 
BR(c -+£) = (9.73 ± 0.32)% 


±0.05 
±0.02 
±0.23 


±0.06 
±0.03 
±0.20 


b — > i modelling 

b — > c — > £ modelling 

c — > £ modelling 


-0.01 
±0.02 
±0.03 


< 0.01 
±0.02 
±0.12 
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(x ( e>) = 0.702 ± 0.008, (x ( e>) = 0.484 ± 0.008 pi], obtaining Ag c c - C = ±0.02%, 
Ag^ = ±0.03% and Ag c cS = ±0.04%, Ag^ = ±0.05%, respectively. 

Lepton identification efficiency 

The charm particles produced from gluon splitting are generally of low energy and 
this has two consequences for the lepton selection. Close to the lower momentum 
cut-offs the identification efficiency for the electrons degrades due to the poorer 
relative energy resolution of the calorimeter, whilst for the muons there is an 
increasing level of backgrounds from ir and K decays, as discussed in (vi). The 
systematic error associated to the lepton identification efficiency is estimated as 
in [IB]. The resulting effect is Ag e c - C = ±0.18% and Ag^ = ±0.01%. 

Fake and non-prompt leptons 

A light hadron (K, it), or a light hadron decaying to a muon within the tracking 
volume (K — > //, tt — > /i), can be selected by the muon identification algorithm 
with a certain probability. A reliable estimate of such a mistag probability is 
important for the evaluation of the purity of the selected sample, derived from 
simulated events. 

Two high-purity samples of light hadrons are selected, and the muon selection 
efficiencies in data and simulation are compared. 

The first sample is selected inclusively by means of a dedicated "uds" tag, 
designed to select event hemispheres that do not contain decay products of heavy 
flavoured particles. Such a tag is based on the presence of secondary vertices, on 
the momentum of the fastest charged particle, on the total visible energy and, if 
a lepton candidate is present, on its transverse momentum with respect to the jet 
axis. Tracks are selected in the opposite hemisphere with the muon identification 
kinematic cuts as in Section El 

A second sample is obtained by reconstructing Kg decays, as in [23] • The 
candidates are required to have a decay length larger than 1 cm, and an invariant 
mass in a window of width (0.014 ± 0.3Pk/\/s) GeV/c 2 around the nominal 
Kg mass. The contamination of prompt muons from heavy flavour decays is 
further reduced by applying a soft uds tag to the hemisphere opposite to the Kg 
candidate. 

The number of tracks selected and identified as muon in data and simulation are 
shown, for the two samples, in Tables El and The mistag probability in the 
data is calculated after subtracting the prompt muon component estimated from 
the simulation. 

The two samples consistently indicate that the mistag probability is higher in data 
than in the simulation, by the ratios 1.17 ± 0.03 and 1.15 ± 0.09. A correction 
factor of 1.17 ± 0.05 has been conservatively applied to the simulation. The 
resulting uncertainty on the measurement is Ag^ = ±0.50%. 

For electrons an uncertainty of 20% is assigned to the rate of misidentified 
hadrons, and of 5% to the rate of photon conversions ^Hj, yielding 
Ag^ = ±0.09%. 
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Table 6: Performance of the prompt-muon identification on uds-fragmentation hadrons. 





Simulation 


Data 


Number of tracks 


564396 


332103 


Number of prompt /i 


432 




Tracks identified as \x 


3072 


2077 


Prompt /i identified as fi 


351 




Mistag (10~ 3 ) 


4.83 ±0.08 


5.64 ±0.13 



Table 7: Performance of the prompt-muon identification on pions from Kg decays. 





Simulation 


Data 


Number of tracks 
Number of prompt \l 


113068 
199 


50162 


Tracks identified as \x 
Prompt fi identified as \x 


662 
156 


326 


Mistag (10^ 3 ) 


4.48 ± 0.20 


5.13 ±0.31 



vii. Generators 

The modelling of the g — > cc process affects the selection efficiencies and hence the 
value of g C c- The process can be described in terms of three basic variables: the 
energy E g of the gluon, its effective mass m* and the decay angle 9* of the c quark, 
measured in the gluon rest frame. The distributions of these variables as given by 
the JETSET generator are reweighted to match the prediction of the HERWIG 
generator [21], and the difference observed in the measured value of g c5 is taken as 
an estimate of the systematic uncertainty. The procedure yields Ag^ = ±0.05% 
and Ag^ = +0.26% for the effect of E g and m* (which are strongly correlated), 
Ag^ = —0.06% and Ag^ s = —0.12% for the angular distribution. 

viii. Mass of the charm quark 

The mass of the charm quark is taken to be 1.2±0.2 GeV/c 2 . A shift in the charm 
mass results in a variation of the E g and m* distributions. The corresponding 
uncertainties on the measured values are Ag^ = ±0.19% and Ag^ = ±0.24%. 



6 Consistency checks 

6.1 Stability of the neural network cuts 

The stability of the results is checked against the cuts on N e and N^. Figure El shows, for 
the electron and muon subsamples, the variation of the results versus the cuts, together 
with the statistically uncorrelated error. No significant dependence is observed. 
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Table 8: Summary of the statistical and systematic uncertainties. 



source or error 


A 9cc {70) 


A 9cc (/o) 


Statistical error 


0.28 


0.38 


Statistics of the simulation 


0.13 


0.21 


9bh 


0.09 


0.10 


Thrust cut efficiency 


0.07 


0.11 


BR(b -> i) 


0.05 


0.06 


BR(b -> c -> £) 


0.02 


0.03 


BR(c -»■ £) 


0.23 


0.20 


b — > £ spectrum 


0.01 


< 0.01 


b — > c — > £ spectrum 


0.02 


0.02 


c — > £ spectrum 


0.03 


0.12 


Lifetimes of b hadrons 


0.03 


0.03 


(4 C) ) 


0.04 


0.05 


(4 b) ) 


0.02 


0.03 


Lepton identification 


0.18 


0.01 


Lepton background 


0.09 


0.50 


Generators, 9* 


0.06 


0.12 


Generators, E g ,m* 


0.05 


0.26 


m c 


0.19 


0.24 


Total systematic error 


0.42 


0.72 



6.2 Analysis of simulated events 

In the simulation, without reweighting, the fraction of hadronic events with a gluon 
splitting to a cc pair is 1.79%. The analysis applied to simulated events yields 
g% = (1.70 ± 0.21)% and g& = (1.75 ± 0.29)%, consistent with the input value. 

6.3 Shape of N e and 

The shape of the neural network output distributions for the electron and the muon 
samples is shown in Fig. |3J The agreement between data and simulation over 
the whole range confirms that the excess observed in the data is compatible with 
originating from gluon splitting to cc. A fit of g cE to this shape would have therefore 
been feasible and would have (marginally) improved the statistical accuracy of the 
measurement. However, because this shape is expected to be sensitive to details of the 
fragmentation simulation, the already dominant systematic uncertainty would have 
increased accordingly, thus reducing the significance of the measurement. 
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Figure 5: Ag c5 versus the N e , cut. The error bars indicate the uncorrelated 
uncertainties. 

7 Results and conclusions 

The rate of gluon splitting to a cc pair in hadronic Z decays has been measured from 
samples containing an electron or a muon candidate: 

g% = (3.32 ± 0.28(stat) ± 0.42(syst)) % , 
g& = (2.99 ± 0.38(stat) ± 0.72(syst)) % . 

The two results are combined with the correlations between systematic uncertainties 
taken into account using the BLUE technique . The combined final ALEPH result 
is 

g c - c = (3.26 ± 0.23(stat) ± 0.42(syst)) %. 

This result is in agreement with the other measurements [3 El El and with the 
present world average g cZ = (2.96 ± 0.38)% |T9"| . 
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